Full Counting Statistics of Non-Commuting Variables: the Case of Spin Counts 
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We discuss the Full Counting Statistics of non-commuting variables with the measurement of 
successive spin counts in non-collinear directions taken as an example. We show that owing to 
an irreducible detector back-action, the PCS in this case may be sensitive to the dynamics of the 
detectors, and may differ from the predictions obtained with using a naive version of the Projection 
Postulate. We present here a general model of detector dynamics and path-integral approach to the 
evaluation of PCS. We concentrate further on a simple "diffusive" model of the detector dynamics 
where the PCS can be evaluated with transfer-matrix method. The resulting probability distribution 
of spin counts is characterized by anomalously large higher cumulants and substantially deviates 
from Gaussian Statistics. 
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I. INTRODUCTION 

In the past years, there has been a growing interest 
in noise in mesoscopic systems^. Normally, noise is an 
unwanted feature, and, according to classical physics, 
in principle can be made arbitrarily small by lowering 
the temperature; according to quantum physics, however, 
noise is uneliminable due to the intrinsic randomness of 
elementary processes. Furthermore, noise, rather than 
being a hindrance, contains valuable information which 
adds to the one carried by the mean value of the quan- 
tity observed. Simple probability distributions, like e.g. 
the gaussian ones, are determined by the mean values 
and noise. Even though gaussian distributions are ubiq- 
uitous, there are interesting physical processes which are 
described by non-gaussian distributions. Noise alone is 
not sufficient for the determination of such distributions. 
One needs to know all the momenta, or equivalently their 
generating function. Full Counting Statistics^''^ consists 
in determining the latter. 

The FCS approach has been receiving increasing atten- 
tion from the physics community. Its connection with the 
formalism of non-equilibrium Green functions^' and cir- 
cuit theory^ was established^. It has been used to char- 
acterize transport in heterostructures^, shuttling mecha- 
nismSiS, charge pumpingiS., and multiple Andreev reflec- 
tionsiiiiS. The technique was extended to charge counts 
in multiterminal structuresi^, and to spin countsi^. The 
FCS of a general quantum variable was studied, and 
the necessity of including the dynamics of detectors 
stressedi^. 

There are some open issues in FCS. The main one con- 
cerns whether it is always possible to find a generating 
function which allows an interpretation in terms of prob- 
abilities. Indeed, in Ref. it was found that such an 
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interpretation is not straightforward. This problem was 
shown to amount to the long standing question of the 
non-positivity of the Wigner distribution^-^'. The lack of 
a classical interpretation was attributed to the breaking 
of gauge invariance for the charge degrees of freedom, due 
to the presence of superconducting terminals.-^' It is inter- 
esting to consider a more general mechanism of gauge in- 
variance breaking which involves spin degrees of freedom. 
It may be caused either by the presence of ferromagnetic 
terminals or by subsequent detectors measuring different 
components of the spin. We shall consider the latter case. 

Another issue we want to address is the range of ap- 
plicability of the Projection Postulate. Since von Neu- 
mann's classic work^^, it is known that Schrodinger's 
evolution cannot account for the fact that the result of 
an individual measurement has a unique value, and can- 
not be described by a superposition. It is necessary to 
supplement Schrodinger's evolution with an additional 
evolution (type II in the terminology of von Neumann) , 
projecting the state of the observed system into the eigen- 
state of the measured observable corresponding to the 
actual outcome. This can be done at several stages: one 
could dispense with the description of the measurement, 
and project the wave-function of the system. Alterna- 
tively, one may continue the chain by describing the in- 
teraction of system and detector, trace out the system's 
degrees of freedom, and then project the state of the 
detector. This chain can be continued indefinetely, by 
skipping the projection of the detector's state, and con- 
sidering the coupling of the detector with the visible ra- 
diation, of the latter with the eye of the observer, etc. So 
far, it has been implicitly assumed that the predictions 
of Quantum Mechanics do not depend on the stage at 
which one chooses to stop the chain and project. In this 
work, we shall demonstrate that different statistics are 
predicted when one projects at the level of the system 
and at the level of the detector. The reason for this is 
that, in the example we shall discuss, the quantum dy- 
namics of the detectors cannot be neglected, even after 
accounting for decoherence. 
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Besides being the simplest illustration of noncommut- 
ing variables, detection of spin components is a worthy 
subject in its own right. Spintronics, i.e. the study of 
how producing, detecting, and manipulating spins, is a 
rapidly growing field^^, which has already found impor- 
tant technological applicationsi^. 



In this paper, the subsequent detection of non- 
commuting variables is discussed. The Full Counting 
Statistics approach allows to obtain the joint probabil- 
ity distribution for the counts. The non-commutativity 
of the observed variables manifests itself in the fact that 
the back-action of detectors, and their quantum dynam- 
ics, must be taken into account. This remains true also 
when an environment-induced dissipative dynamics for 
detectors is included. The reason is that one does not 
observe one particle at a time, but a flux of particles 
traversing the detectors at a rate which can be larger 
than the decoherence rate of the detectors themselves. 



The present paper is laid out as follows: In section UTI 
the connection of FCS with the density matrix of the de- 
tectors is derived, and a general theory of detection of 
non-commuting variables is presented; a model for the 
measurement is introduced. In section IIIII we discuss 
the case of ideal quantum detectors, having no internal 
dynamics. We argue that they do not provide a realis- 
tic model of detectors because of their long memory. In 
section Hvl we discuss the internal dynamics of the de- 
tectors. The fact that detectors are "classical" objects is 
accounted for by introducing a dissipative dynamics due 
to their interaction with an environment. Then, since 
we intend to concentrate on spin counts, in section Ivl we 
present a model for a spin detector in solid state, relying 
on spin-orbit interaction. We proceed to section IVII by 
introducing the particular system that we study, namely 
two normal reservoirs connected by a coherent conduc- 
tor. In section lVIll we give details about the derivation of 
the FCS for this system, relying on the full quantum me- 
chanical description of detection process, and we present 
the results. In section IVIIII we discuss the FCS that 
would be obtained by a naive application of the projec- 
tion postulate, i.e. by neglecting the quantum dynamics 
of the detectors. In section Hxl we compare the results of 
the two approaches for the case of one and two detectors 
in series, and we find that they coincide. In section ^ 
we find a discrepancy between the two approaches when 
three detectors in series are considered. In particular, we 
show that both approaches predict the same second-order 
cross correlators, and that they differ in the prediction of 
fourth order cumulants ((crjcrl)). Finally, in section IXll 
the case of three spin detectors, monitoring the X, Y and 
Z components of spin current, is presented. The proba- 
bility distribution for the counts reveals large deviations 
from the Gaussian distribution. 



II. GENERAL CONSIDERATIONS ABOUT 
MEASUREMENTS 

All the information that we can gain about a system is 
stored in the density matrix of one or more detectors (de- 
noted by index a) which have interacted with the system 
during a time r. The reduced density matrix is 

where Tr^j^^ stands for the trace over the degrees of free- 
dom of the measured system, p{0) is the initial density 
matrix of system and detectors, and Urfi is the time 
evolution operator. We focus on the representation of 
Pdet in a basis \(p), p'j^'J^ (r) ee {(l)\pdet{T)\4''). Here, 
10) — \(j)a) is a vector in the Hilbert space of the 
detectors. Since the time-evolution is linear, a matrix 

zt't, exists such that 

Ptt'ir)^ j dpdp'Zllip^XQ). 

Thus, given that one knows the initial density matrix of 
the detectors, Z contains all the information one can ex- 
tract from the measurement. However, part of this infor- 
mation gets lost: we can only know the diagonal elements 
of the density matrix in a particular basis, identified by 
the pointer states of the detectors. These states, which 
will be denoted by |7V), correspond to the detectors in- 
dicating the values {Ng}, and are individuated by the 
property that, if one prepares the detector in a generic 
state identified by a density matrix p^^^ , and then lets 
the environment act on it, the off-diagonal elements of 
the density matrix in the basis |A'^) will go to zero with 
an exponential decay. We point out that this does not 
dispense us from invoking a projection at some point. 
The presence of the environment explains how the en- 
semble averaged density matrix reduces to diagonal form 
in the basis of pointer states, but it does not explain 
how the density matrix of the subensemble correspond- 
ing to an outcome Na purifies to the state \Na). This 
requires invoking the projection postulate for the detec- 
tor, or, equivalently, an evolution dictated by the rules of 
the bayesian approach^^ or of the quantum trajectorySS 
one. 

The quantity accessible to observation is the probabil- 
ity to find the detectors in states \Na), after a time t. It 
is given by 

PAN) - {N\paet{T)\N) . (1) 

If offdiagonal elements of the detector's density matrix 
decay instantaneously, Pt {N) depends only on the prob- 
abilities at a time immediately preceding r, Pr-dtiN), 
and the process is Markovian. 

Now, let us consider the operators Ka corresponding to 
the read-out variables of the detectors. Their eigenstates 
are \Na), where Na indicates an integer which is propor- 
tional to Ka- The proportionality constant is provided 
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below. Let us also introduce the conjugated operators 



Va 



Ka,Vb 



iSabfi, and their eigenstates 



, with 

dimensionless quantitities proportional to Va- If we in- 
sert to the left and to the right of pdet in the RHS of 
Eq. the identity (in the detectors' Hilbert space) in 
the form I oc J we obtain 



2n 2n 



X exp 



(2) 



We used the shorthand (f) ■ N = J^a 't'a^a- We change 
variables according to 0=*= = (<i> ± (/>)/2. Here, $ and (f) 
are the classical and quantum part of the field, respec- 
tively. This terminology reflects the fact that fluctuations 
of $ are set by the temperature, while fluctuations of (j) 
depend on S, as we shall prove in section Hvl The time 
evolution depends on the Hamiltonians of the system and 
the detectors, and on their interaction. We focus on the 
detection of internal degrees of freedom of a system whose 
center of mass coordinate x is not affected by the pres- 
ence of the detectors. We consider several detectors in 
series along the path x(t) . We take the interaction to be 
of the form Hi^t = Y.a ^int with 



Htnt = -aa{^)XaVaJa , 



(3) 



where x is the coordinate of the wave-packet, Xa cou- 
pling constants depending on the actual detection setup, 
Ja is an operator on the Hilbert space of the system's de- 
grees of freedom, and aai'x.) is a function which is unity 
inside the sensible area of the a-th detector and zero out- 
side. For a one-dimensional motion, e.g., we would have 
a{x) = 0{x-X'^"'y)e{X'^f"'^ -x), with and X'^f"'^ 

are the coordinates delimiting the sensible area of the de- 
tector, 9{x) = if a; < 0, d{x) = 1 if x > 0. Ja is the 
current associated with the measured quantity, such that 
the output of the detector does not depend on the time 
each particle takes to cross its sensible area. Indeed, the 
equation of motion for the "measuring" operator is 



dKgjt) 

dt 



aa(x(i))AJa(t) 



In the equation above, we have assumed that the operator 
Ka commutes with the unperturbed Hamiltonian of the 
detector. In general, however, < Ka > will fluctuate 
in time due to background noise. Such fluctuations put 
a lower limit to the resolution of the detector. For a 
reliable detection, the resolution must be smaller than 
the minimal variation Kqa one intends to measure. 

Let us introduce proper units. We consider the case 
where the measured quantities have discrete values pro- 
portional to a quantum Eqa- For instance, for charge 
Eq = e, the elementary charge, and for spin Eq = h/2. 
Every time an elementary unit passes the detector, the 
readout of the latter will change by Kga = XaEQa- 



Thus, we introduce the number and phase operators 
Na = KjKQa, = V/VQa, with Vq^ = h/ Kq^. We 
further assume that i) the detectors are initially prepared 
in a state with zero counts pdet{^) = |iV = 0) (iV = |, 
and ii) the spread of the system wave-packet is much 
smaller than the distance between two subsequent detec- 
tors, Ax ^ Xa+i — Xa- The first assumption implies 
that 



ptt(r)=Z(0,0') 



2tt 2tt 



or, explicitly, 



dpi~^ dp 

~2j^^2nx' 

fj.+ 11- 

exp (S'det - Sdet[<P^] + Tsys[(t)^ A']), 

(4) 

where the limits of the path-integrals fix the values of the 
fields a.tt — Q and t — t, and we introduced the influence 
functional of the system on the detectorgSi 



sys 



exp [Tsys[4) \ 4) ]) := Tr 
|exp (^Sint[4>^ ,J'\jpsys{^)ey^T> (-'S'mt[0", J]) |, (5) 

where Sint is the action corresponding to the interac- 
tion Hint given in Eq. (PJ. We shall call Z the quantum 
generating function. In principle it depends on twice as 
many parameters than the classical generating function 
does. In the rest of the paper we shall use the cumulant 
generating function (CGF), T = logZ . The advantage 
of working with the CGF is that it often has a clearer 
interpretation than , since independent processes con- 
tribute factors to Pr and simply additive terms to the 
CGF. Hence, if subsequent events are independent, the 
CGF is proportional to the observation time r. Thus, 
time averaged cumulants, which for long r correspond 
to zero-frequency noise and higher order correlators for 
currents, have a finite value. 



III. DETECTORS WITH NO DYNAMICS 

We analyze the situation where the dynamics of the 
detectors is neglected. This means that 



i.e. the counting fields are constant. 

^e consider first the case of one detector. Then 



^Trsyslu'^^Psys{0)U'^~^}, (6) 
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where 



l(f> / dtJ{t)/EQ 



(T being the time-ordering operator) is an operator in the 
system's Hilbert space. By exploiting the cycHc property 
of the trace, we have that, if J is a conserved operator or, 
more generally, [ J(t), J{t')] = 0, then Z{(l)~^ , (f>~) depends 
only on (j) — 0+ — 0^ . It has been shown that in this case 
Z(0^,0+) gives directly the generating functioni^. 

Next, we consider the case of two detectors. The kernel 
2 is now 



Tr. 



(7) 



Here we exploited assumption ii), and defined 



exp 



l^a / dtJa{t)/EQ 



Once again we exploit the cyclic property of the trace and 
see that the expression does not depend on the combina- 
tion $2 = '/'2' + "^2^ • From Eq. ^ we see that in general, 
for two detectors, Z does depend on <I>i, even when Ja 
are conserved. However, when the system is initially in 
the unpolarized state psys oc Isys, the dependence on $i 
disappears as well. Another case in which this happens 
is when the detectors monitor two commuting degrees of 
freedom which are conserved. For instance, if the current 

J is not conserved, in general J{t),J{t') ^ 0. Thus, 

even if one repeats the same measurement, one would ob- 
tain different results. If however the current is conserved 
J{t) — J, and both detectors measure J, the kernel de- 
pends only on the combination (j)i+(j)2, which means that 
the two measurements will give the same result. 

In general, when there are three detectors, labelled 
1, 2 and 3 according to their order, measuring non- 
commuting quantities, even if the system is initially un- 
polarized, the integrand will depend on the classical vari- 
able of the middle detector, $2- When such a dependence 
appears in the expression for the generating function, it is 
a signal that the internal dynamics of the detector must 
be taken into account. Indeed, when Z does not depend 
on $, the density matrix is diagonal in the basis |A^). 
When Z does depend on <&, pdet develops off-diagonal 
components. We consider as an example the case where 



the detectors' density matrix is prepared in a diagonal 
state at i = 0, and two particles are sent to the detec- 
tors one at time ti > and the other at time t2 > ti, 
in such a way that their wave-packets do not overlap. 
Then, after the first particle has crossed the detectors, 
the density matrix of the detectors p^^^ has offdiagonal 
elements, which depend on the original diagonal elements 
(probabilities). Since one observes only the probabilities, 
this can not be ascertained directly. However, when the 
second particle crosses the detectors, the new probabil- 
ities will be a combination of the former diagonal and 
off-diagonal elements. In order to know P^i^^ {^2), knowl- 
edge of p^^^ (ti) is not sufficient. Thus, the process is 
non-Mar kovian. In principle, even after the detector has 
been measuring for a long time a large number of parti- 
cles, the off-diagonal elements created after the passage 
of the first particle will still influence its dynamics. This 
is not realistic, since, because of the coupfing of the de- 
tectors to the environment, the offdiagonal elements will 
go to zero within a typical time Tc- In order to account for 
this, one should consider the dynamics of the detectors, 
which we shall do in the next section. 



IV. DETECTORS WITH INTERNAL 
DYNAMICS 

We model the decoherence of the detectors by intro- 
ducing a dissipative dynamics for the detectors' degrees 
of freedom, i.e. we couple the detectors to an environ- 
ment, whose degrees of freedom are traced out. We model 
the environment as a system of independent harmonic os- 
cillator in thermal equilibrium, having the action 



2 2 21 



and coupling to the detectors through the position oper- 
ator 



Sdet-env — Jl ^ ^ 



CjaX]VQa4>a, 



(8) 



with Cja coupling constant between the j-th oscillator 
and the a-th detector. Then the generating function be- 
comes 



0+ Xj Xj 

Z{4>+,r) = j / dxjdx+dx- Jvc^+{t)Jvr{t) jvx+{t) Jvxj{t)p,Ux+,x-) 

M- 

expjs'deti^"'"] - Sdet[4>~] + Fsys[(t)~^ , + Senv[x^] - Senv[x~] + Sdet-env[x^ , (j)^] - Sdet-env[x~ , , (9) 
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In the expression above, we isolate the part 

Xj Xj 

expj^env — J dxjdx'j'dxj J'Dxj'{t) Jvxjlt) 

x'^ 

J J 
- Sdet-env[xJ , (1)^]^ Penv{x^ , X^) (10) 

which gives the influence functional of the environment 
on the detectors. We notice from Eq. ||SJ| that, since the 
functions (/>a (^) fixed by the external path-integrals, 
they act as an external source ifit) = J2a'^ja.4'ti^) 
the j-th harmonic oscillator. It is then possible to per- 
form the independent gaussian path- integrals over Xj , re- 
sulting IdSS 

^e™ = - ^ E ^Qa [ dt dt' - 0- (t)) 

x[aa{t~t')^Ut')'C.:{t^t')^-{t')] , (11) 

where the influence of the environment is contained in 
the complex functions aa{t), whose Fourier transforms 
are 

aa{uj) - i (^coth + 1) <ya{uj) , (12) 

where the inverse temperature /3 = l/ksT comes from 
having assumed the bath in thermal equilibrium {penv — 
exp {—f3Henv)), and CTq are the spectral densities 



At low frequencies, we can approximate the odd- 
functions (7a by (7a (w) ~ (Ohmic approximation), 
with 7a friction constant, as will be clear later. We in- 
troduce new variables (j> — 4>~^ — , $ = -I- (/)~ . Thus 
we get 

1 f dcu Bhuj , Bhiu , , , ^,1 
-^J^'—^ot^^\U^)f (14) 

We take the action of free detectors to be that of har- 
monic oscillators, i.e. 



S, 



det\ 



where uia is the "mass" of the detector (i.e. it is 
the inertial term corresponding to the kinetic energy 



moVg^^a/^)- Then the generating function reads 



M 



exp 



E 



2h 



Zn 



la 



where we introduced the response function 



2 , ■ 7a 

m„ 



(16) 



(17) 



'a 



from which one can see that 7a are proportional to the 
friction constant, and the fluctuation term 



— — — coth— — . 



(18) 



The part of the action containing the fluctuation term in 
(j){uj) is, at low frequencies, proportional to temperature 
T and to l/ti^. The factor l/fi^ strongly suppresses large 
fluctuations in 0. Thus, the influence functional due to 
the measured system J^sys[^,(l)] = J dtLinf{^{t),(j){t)) 
can be approximated by J dtLinf{^{t),(t)) := ^<t>[^]- In- 
tegration over <pc^ gives finally 



Z(0,$) =y dM jvm) e^*[*l+^=//[*l 



(19) 



M 



with the effective action 



2 \2 



1 ^ {l3maV^^ 



la 



duj 



\9a\u;W 



\^aiL0)\' 



We notice that at high temperatures /(i^) — 1, and thus 
h disappears in the effective action for <&. For this reason 
the latter is termed the "classical" part of the field. 

In the limit of small mass TOq 0, niai^^VQ^ Ea, 
where Ea has a finite value and is a typical energy scale 
of detector a, the effective action simplifies to 



Seffm 



dt 



^a{t) 



1 



with Tac = P'yaVga/^ the "coherence time" of the detec- 
tor, and A$a = 2/PEa the spread of $a. 

V. SPIN DETECTOR 

We discuss a model for spin detection. The setup cor- 
responds to the one proposed and used in~i to detect 
Aharonov-Casher effect for neutrons. This setup ex- 
ploits the fact that a moving magnetic dipole generates 



6 




FIG. 1: The proposed spin current detector. An electron 
with velocity v and spin S induces a voltage drop in a ca- 
pacitor. The electric field E inside the capacitor produces an 
Aharonov-Casher phase shift on the electrons. 



an electric on o^^i^^ . To measure this, one encloses the 
two-dimensional current lead between the plates of a ca- 
pacitor as shown in Fig. ^ While in Ref. the neu- 
trons passed a fixed electric field, which gave a constant 
Aharonov-Casher phase, in a spin detector the initial 
voltage applied to the plates is zero, and the passing of a 
particle with spin 1/2 will cause the charge in the capac- 
itor to show pulses towards positive or negative values 
depending on the result of the measurement. The associ- 
ated phase Kt = /J dtQ{t) will thus increase or decrease 
stepwise in the ideal situation where spins are transmit- 
ted separately in vacuum through the detector. 

Each spin moving with velocity v produces an elec- 
tric field. For electrons in vacuum, the interaction term 
between spin and detector is given by the spin-orbit cou- 
pling 



1 ^ , V . 



with c the speed of light, and the factor 1/2 accounts 
for the Thomas precession. The magnetic moment fi 
is proportional to the spin fj, = {gs\e\/2me)S, with nie 
mass of the electron, e = — |e| its charge, and gs its spin 
gyromagnetic factor. Thus, we rewrite the interaction as 

Hm - -{gs\e\/4.m,c^)E ■ {v x S) . 

The spin-orbit coupling induces a current in the RC cir- 
cuit. The integrated charge traversing the circuit is the 
detector read-out. The read-out signal is proportional 
to spin current in the lead 3, Q = An • J, n being the 
unit vector perpendicular to the direction of the current 
flow and parallel to the plates of the capacitor, A being 
a proportionality coefficient. The concrete expression for 
the latter, A = gs\e\L„/ AuieC^w, depends on the geomet- 
rical dimensions of the detector the length of its plates in 
the direction of the current L,,, and the distance between 
the plates w. The variable canonically conjugated to the 
read-out is the voltage V across the capacitor, and the 



expression for the interaction in terms of V contains the 
same proportionality coefficient A, Hint = —XVn-3. Our 
choice of the detection setup is motivated by the fact that 
this detector does not influence electron transfers through 
the contact and only gives the minimal feedback compat- 
ible with the uncertainty principle: the electrons pass- 
ing the capacitor in the direction of current acquire an 
Aharonov-Casher phase shift, which consists in a preces- 
sion of the spin around the detecton axis n. This depends 
on spin and is given by ^ac = XVn-S/h. This is similar 
to the detection scheme presented inm for charges trans- 
ferred. A fundamental complication in comparison with 
the charge FCS is that in our case the phase shift depends 
on spin, so that even the minimal feedback influences 
the statistics of the outcomes of following spin detectors. 
We introduce dimensionless variables N = 2 J dtQ/hX, 
(f) = XV/2. Then N varies by one every time a spin 1/2 
crosses the detector. With reference to Eas. (|17|) . (|18f) . we 
have m = LC^ 0, = 1/LC oo, mfl^ C, 
7 RC'^, E = 4C/A^, with L,, R and C inductance 
(assumed negligible), resistance and capacitance of the 
circuit. 



VI. THE SETUP CONSIDERED 

We consider a system composed of two metallic, un- 
polarized leads, connected through a coherent conduc- 
tor, characterized by a set of transmission probabilities 
r„, where n identifies transmission channels. A negative 
bias voltage V is applied to the left lead. At the right of 
the conductor there are several spin detectors, labelled 
from left to right by a = 1, 2 • • • , and a current detec- 
tor, denoted by a = 0. The counting fields will be then 
(j)a, with a = 0, 1, • • • . Since charge and spin currents 
commute, the current detector can be positioned at any 
point along the chain of detectors, without influencing 
the statistics of the outcomes. The setup is depicted in 
Fig. |5| We require that the coherent conductor is non- 



Contact 
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FIG. 2: The setup considered, in the case of three spin detec- 
tors and one charge detector. 

polarizing. Thus, the average spin current is zero. How- 
ever, there are spin fluctuations, which are revealed by 
measuring noise and higher order correlators (or cumu- 
lants). From the symmetry with respect to reversal of 
spin, we can predict a priori that all odd cumulants are 
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zero. We shall concentrate on a situation where there are 
three spin detectors. This is because, as anticipated in 
section ^1 the current is unpolarized and one needs at 
least three detectors monitoring non-commuting quanti- 
ties in order to see non-trivial consequences of the detec- 
tors' feedback on the system. The feedback consists in 
the wave-function picking up an Aharonov-Casher phase 
while traversing each detector. 



Gn, we obtain 
Gr^ 

where 



1 - -2fRe'^M 
-2(1 - ,fe)e-'x^t 2fR - 1 



^-11' 



-i(0-/2)n„ 



(21) 



VII. RESULTS 

The technique we use is an extension of the scattering 
theory of charge FCS. This theorySii^ expresses FCS in 
terms of a phase factor e*-^ acquired by scattering waves 
upon traversing the charge detector. 

Since we do not consider energy-resolved measure- 
ments, the phase factor does not depend on the channel, 
and the approach works for a multi-channel conductor 
as well as for a single-channel one. The phase factor e*-^ 
can be seen as resulting from a gauge transform, to be 
applied to the (known) Green function of the right lead, 
that removes the coupling termSi^ Hint = ~f ■ For the 
case of the spin detectors, the gauge transform introduces 
a phase factor which is a unitary matrix in spin space. 
Namely, the gauge transform generated by spin detector 
a is e*'^""" '^. In this matrix, r is a pseudovector of 2 x 2 
Pauli matrices, and Wa is the direction along which detec- 
tor a detects spin current. The Keldysh Green function 
of the lead is 



Gi{E) 



1 - 2/, -2fi 
-2(1-/0 2/,-l 



where I E {L, R} denotes the left or right lead, and // is 
the corresponding Fermi occupation number at energy E 
and chemical potential The elements of the matrix 
are actually in their turn a matrix in spin space. Since 
the leads are assumed to be unpolarized, they are simply 
the identity. The matrix current is given bjii 



T 


Gl , Gr 




l+Tn 


({G'L,Gfl,}-2) 


/4 



(20) 



from which it follows that the quantum generating func- 
tion is 

^-^E/^^log(l + ^"({6.,G4-2)/4). 

n •' 

Here [...] ({...}) denote (anti)commutator of two matrices, 
and Gr is the transformed matrix 



Gr — 



rR = e-- I I e 

a 



-Upalla-T ~IX 



where x = diag(x"'", x ), = diag (0+, ip ) are matrices 
in Keldysh space. After substituting the expression for 



is a matrix in spin space. We notice that i) the charge 
fields come only in the combination x = x+ — x~ , and 
ii) the phase factors e*-^ in the expression for charge FCS 
are replaced by e^^Ai to give the FCS of charge and spin 
counts after taking trace over spin. If we also notice that 
is a (2 X 2) matrix with eigenvalues e^'", we arrive to 



•^(X,{0^},{0a}) = JE-^-(^±")' (22) 



where ^dx) is the generating function for charge count- 
ing. The a is given by 



cos a = ^trM 



(23) 



The explicit expression for the system considered here, 
in terms of the transmission probabilities through the 
contact and the applied bias is, at zero temperature, 

^ = /■ ^ y log [i?2 + T^^x + 2RnTne'^ cos a] , 

with Rn — 1 — Tn, TV = 2Trh/eV. The interpretation is 
quite straightforward: electrons coming through differ- 
ent channels behave independently, which is revealed by 
the fact that the generating function splits into a sum; 
each channel can accommodate two electrons in a spin- 
singlet configuration; with probability R^ none of the 
two electrons passes the junction, and there is no contri- 
bution to the charge counting nor to the spin one; with 
probability T^^ both electrons come through the conduc- 
tor. This gives a contribution of two elementary charges 
transferred (factor e^*-^), but no spin transfer. Finally, 
with probability p„ = 2_R„r„, exactly one of the two 
electrons is transferred. This gives a contribution to the 
charge and to the spin counting. 



VIII. PROJECTION POSTULATE 

We demonstrate that a different FCS is predicted by 
using a different approach, namely a naive application of 
the projection postulate, consisting in avoiding the de- 
scription of the measurement and applying the projection 
to the system measured. We shall denote this procedure 
with PP for brevity. This approach predicts a parameter 
app which does not depend on $. Let us give the details 
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of such a derivation: When an unpolarized electron ar- 
rives to the first detector, the probabihty of the outcome 
CTi = ±1 is Pi((Ti) = 1/2. The conditional probability 
that the second detector gives (72, given that the first 
read cti is P2(o'2|o'i) = (1 + (Ti(T2ni • n2)/2. This is be- 
cause after the first detection the spin of the electron is 
assumed to have collapsed along ±ni. The same happens 
after the second detection. Consequently, the conditional 
probability that a third detector reads crs, given that the 
first read tri and the second <ji , depends only on the lat- 
ter outcome P3(o-3|o-2, cti) = (1 + cr2CT3n2 • n3)/2. The 
process is in a sense a Markovian one. The total joint 
probability for each electron transfer with an arbitrarily 
long chain of detectors is 

P(W)=2ll ^ , 

a=l 

and the corresponding generating function for the setup 
considered here is given by Eq. H24|l with 



the read-out of the charge is not exploited (x = 0) . Then 



cos app = cos 



K-l 



\-\- (TaCFa+ina-ria+l 



(25) 



IX. COMPARISON OF THE TWO 
APPROACHES FOR ONE AND TWO SPIN 
DETECTORS 

Now, let us go back to Eqs. if^ . lt^ and compare 
the two approaches for some simple cases. For the case 
of one or two detectors in series, the eigenvalues e"*"*" 
are not affected by the order of matrix multiplication in 
(|21|l and depend on differences of spin counting fields 
0a = 0+ — 0a only (in fact they coincide with the value 
giiQpp) This impUes that the FCS definition (g)) can be 
readily interpreted in classical terms: it is a generating 
function for probability distribution of a certain number 
of spin counts ctq in each detector. 



^(K}) = j Wdcjya 



(26) 



For a single detector, the spin FCS is very simple: it 
corresponds to independent transfers of two sorts of elec- 
trons, with spins "up" and "down" with respect to the 
quantization axis. The cumulants of the spin (charge) 
transferred are given by the derivatives of F with respect 
to 01 (x), at X = 01 = 0. In this case a = 0i. From this 
and relation (|22|l . we conclude that all odd cumulants of 
spin current are 0, as anticipated, and all even cumulants 
coincide with the charge cumulants. 

For two spin detectors, with rii • n2 = cos0, we obtain 
cos a = cos 01 cos 02 — sin 0i sin 02 cos 9. Since there is no 
dependence on <i>Q, the quantum generating function has 
an immediate interpretation; we consider the case when 



-2(0) = J]^ [q-n + Pn COS a]^\ 



(27) 



where Pn ~ 2i?„T„ is the probability that, in two at- 
tempts of transmitting one electron over a spin degener- 
ate channel n, exactly one is transmitted and 5„ = 1— Pn- 
This result coincides with what one would obtain from 
the Projection Postulate. 

We discuss in detail the probability distribution. By 
performing the Fourier transform, we find the probability 
of detecting a spin ui in direction ni and 02 in direction 
n2: 



p(cTi,a2)=5]'^'n^«(^i"',^ 



-^")) , (28) 



where the prime in the sum means that it is restricted to 
X^n '^^'^ ~ ' ^'^'^ probability for each channel n is 

Pn{cTi,a2) = Ptr{k\N)Py{{k + (7i)/2|fc) 

k 

X I]A((fc + ^2 + 20/4|(fc + ai)/2 =r) 



xP|((A: + (72 - 20/4|(fc-(7i)/2 =i) 



where 



Ptr{k\N) = ( ]plq^~'^ 



N 



1 fk 



(29) 

(30) 
(31) 



pm =T) = ( ) [cos^ {omw iomf-' , (32) 

P^{l\k=i)= ( J)[cos^(0/2)] 



lfc-'r„;„2 in /oMi 



(0/2)]' (33) 



The sums are over all values for which the binomials 
make sense (no negative nor half-integer values). Thus 
k,l,ai,a2 have the same parity. P can be interpreted as 
follows: since the current is unpolarized, we can think 
of it as carried by pairs of electrons in singlet configu- 
ration. Then, there is a successful attempt to transfer 
spin when exactly one of the two electrons is transmit- 
ted. This gives Ptr{k\N), the probability of transferring 
k spins over N attempts (p„ probability of success for a 
single attempt) through channel n; the second binomial 
comes from the ways one can pick A^ij — {k + ai)/2 spins 
up out of k spins, with probability 1/2 (we recall that 
the incoming electrons are unpolarized); the third term 
comes from the fact that, given that Ni^ — {k + ai)/2 
spins up according to the first detector are passed to the 
second one, the latter will measure (k + a2 + 2Z)/4 of 
these as spins up (the probability of agreement between 
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detectors being pag = cos^(6'/2)), and the rest as spins 
down; analogously, the latter term comes from the fact 
that given Nii = (fc — (Ti)/2 spins down along direction 
Ml have been detected, {k + (T2 — 2Z)/4 of them will be 
detected from the second detector as spins up, while the 
remaining ones will be detected as down. 

When the two detectors have parallel orientation {9 — 
0) , the second sum in Eq. (|29|l is nonzero only if ai =02, 
giving 

k 

i.e. there is perfect correlation, as is to be expected. 
When the two detectors have orthogonal orientation [6 = 
7r/2), it is possible to perform analytically the sum over 
m: 

P(ai, a2) = P{k\N)P^{{k + ai) /2\k)P^{{k + cj2)/2\k) 

k 

i.e. the outcomes are independent, given that k successful 
spin transfers happened. 



X. COMPARISON OF THE TWO 
APPROACHES FOR THREE SPIN DETECTORS 



where C = cos6'i2 cos023, and the first term is the PP 
result. The second term, as expected, has a typical signa- 
ture of interference effects: it is suppressed exponentially 
if the variance of the corresponding Aharonov-Casher 
phase ($2) ^ 1- Since ^ac is inversely proportional 
to ft, this is the classical limit. In this limit, the result 
coincides with the PP. 

However, fourth cuniulants show a large deviations 
from the PP result. Namely: 

{{<yl<yl)) = {{<yl<yl\p + s^^^iv^))^, (36) 

where A = sir? 612 sin^ 623 , and the PP result is expressed 
in terms of charge cumulants as 

ii^l^l))pp= I [(1 + 2C2)((iV4)) + 2(1 - C'){{N'))] . 

This deviation results from correlations of $2 at time 
scale Tc- To estimate the result, we notice that the charge 
cumulants are of the order of T/Tg,i, t^i being the average 
time between electron transfers. It is easy to fulfill the 
condition Tei Tc <^ t, and in this case {{cf\a'^Jj is much 
larger than PP result. 

It is interesting to study further the probability distri- 
bution which gives rise to such anomalously large fourth- 
order cumulants. This we shall do in the next section. 



For the case of three detectors, we have 

cos a — cos app ~ sin 612 sin 6*23 sin ($2 — ^^2^) sin (^3 sin 4 
cos app = cos 4>i cos (j)2 cos ^3 + 

— cos 9i2 sin (f)i sin 4>2 cos 03 — cos 623 sin 02 sin 03 cos 0i — 
cos 9x2 cos 023 sin 03 sin 0i cos 02 . (34) 

Here Ojk — arccosrij-rifc are the angles between the 
polarizations (n) of detectors j and k, and cos<I'2°'' = 
(ni X n2)-n3/ sin 6*12 sin 6*23, sin*^"^ = (nixn2) • (n2 x 
n3)/ sin 012 sin 023- As before cosapp is the part corre- 
sponding to the Projection Postulate. We notice that 
when two consecutive detectors are parallel or antiparal- 
lel, then app = a. This is because the same measure- 
ment is repeated twice, and thus we fall back to the case 
of two detectors. 

In general, however, cos a depends on <I>2, and thus one 
needs to account for the dynamics of the second detector 
in order to get the probability distribution for the spin 
counts. We recall that the corresponding detector's ac- 



fluctuations of $2. 



, with Tc 



tion is 5'[<I>2] = / dt\ 

coherence time and 

We have calculated the second cumulants or cross- 
correlators: we found that they differ from the ones ob- 
tained by using PP only by small terms. The correlator 
between first and third detector's readings is: 



(('11^3)) = ((iV2)) \c + (cos 013 - C) e-<*')/2 



(35) 



XI. A PARTICULAR CASE 



We discuss for definiteness the case of three detec- 
tors oriented along three orthogonal directions forming 
a right-handed basis. This implies $2"'' = 0. We con- 
centrate on the joint probability distribution for the out- 
comes of the first and the third detector, irrespectively 
of the reading of the second detector. We consider the 
"classical" limit ($2) 00. Then, the generating func- 
tion for the probability P(cri, (73) for counting ui, tis spins 
in the detectors is: 



/l-<S>2,f ( l-T 

d^2,^d^2J / V^2{t)expl / dt 
J$2.i ["'0 

+ — y"ln[g« +PnCOSa(0,$2)] \ 

i TV ^ 



(37) 



where 



cosa(0, $2) = cosapp — sin 03 sin0i cos $2, 

cos app — cos 01 cos 03 . (38) 

We have a path-integral over imaginary time. We ex- 
ploit the quantum mechanical technique and re-express 
the path-integral in terms of amplitudes: 



^(0) = / d'^2,^d^2J ($2/ ; i = zr, 01 t = 0, , 
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Here the counting fields are parameters, and the time 
evolution of the variable $2 is dictated by \^2',t,'P) = 
e-'^(<^')t|$2;O,0), with the Hamiltonian 



Here, q* is the solution to the equation 



1 

2tc 



TV 



dn [qn +Pn cosa(0, $2) 



(39) 

Then, for large values of r, the path-integral can be ap- 
proximated 



(40) 



where EQ{(j)) is the ground state energy of the Hamilto- 
nian. 

The next step is to find an explicit expression for the 
probability. We recall that the probability to have detec- 
tors 1 and 3 measure average spin currents Ii — ci/t, 
h = <j^/t is related to Z{4)) through 



-ir(0i/i+03/3) 



Since we are in the large r limit, we can evaluate the 
integrals in the saddle-point approximation, and obtain 

P{Iuh) cx exp hSo(0*) - iT{cl>lh + 05/3)] , 
where 0* satisfy the saddle point condition 

+ ih =0 



(41a) 
(41b) 



Assuming that the solutions are much smaller than 1, 
(/)* <C 1, we have that the Hamiltonian can be rewritten, 
including terms up to second order in 0, as 



1 d'^ 

2tc d<^l 



1 

2^ 



20103 cos $2) 



(42) 

where we introduced the average time between spin 
transfers, T5 = We recognize the Hamilto- 

nian for the Mathieu equation 

- - ^ + 2gcos(2i;) . 

Thus the ground state energy depends on the lowest 
Mathieu characteristic function ao{q), with the coupling 
strength given hy q = 4(rc/r5)0i03. Namely, 

Eoicj}) = aoiq)/8r, + (0? + cj)l)/2Ts . 

The saddle-point equations H41|l can then be combined 
to give a trascendent equation for q, from which one ex- 
presses 0*, which are purely imaginary, according to 



i(Pz = Ts- 



h - {l3/2)a'o{q*) 

l~a'o{q*)y4 
I, - ih/2)a'o{q*) 
l-a[,(g*)V4 



(43a) 
(43b) 



(44) 



where we introduced dimensionless currents Va — 
y/T~f^Ia. Eqs. are valid in the limit Tsia <- 1, 

i.e. Va < \fr^Ts- 

Finally, we have that the probability distribution is 

logP(/i,/3)cx-ao(<Z*)/8 

-(^i + ^3)'(l+a^,(g*))/[2 + a;,(<Z*)]' 

(45) 

This probability distribution is to be compared with 
the one predicted by applying the PP. The latter is, in 
the same regime TsIa ^ 1, the independent combination 
of two gaussians: 



logPpp(/i, 73) cx -{vl + vl)l2 , 



(46) 



the proportionality constant (t/tc) being the same. 

In the limit Tc <^ ts, we have that Eqs. H45(l and H46|l 
coincide. However, by taking into account that the de- 
tectors have a finite decoherence time Tc, and that the 
time between spin transfers can be much smaller than 
Tc, we find that the probability distribution deviates sen- 
sibly from Eq. H4()l) . This deviation is larger in the regime 
1 <C — 1 1^3 1 <C ^Tc/ts- when both (dimensionless) 
currents are comparable in module and large with respect 
to 1. 

When z^i ^ 1, we find that 

logP (X -1/2/2 + /(^^sM) , (47) 

with the scaling function f{x) defined by 

ao(go(a:)) 1 
/W = ^ + -^xqoix) 



where the condition 



dap 
dq 



= 2x defines qo{x). In 

particular, f(x) diverges at x = 1 according to f{x) ~ 
— 1/16(1 — 2;). In Fig. Owe draw the logarithm of proba- 
bility as a function of 1^3/ I'l for several values of vi, and 
compare with the probability predicted by making use of 
PP. 



XII. CONCLUSIONS 

We have discussed the Full Counting Statistics of non- 
commuting variables. As a concrete example, we fo- 
cused on spin counts in a two terminal device with non- 
ferromagnetic leads connected through a non-polarizing 
coherent conductor. We have provided a formula con- 
necting the FCS of spins to the one of charge. We have 
seen that it is crucial to have a coherent conductor with 
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FIG. 3; The log of probability as a function of Vi/v\ for 
different values of v\ for the configuration studied in the text. 
All the curves have been shifted by V\/2. The upper curves 
correspond to the result of the PCS approach, and the lower 
ones to the PP. The black dotted curve is the limiting scaling 
curve discussed in the text. 

finite transparency connecting the two leads. This is be- 
cause electrons transmitted through the same channel are 
in a spin singlet, and thus contribute no net spin transfer 
nor spin fluctuations. However, if the transmission prob- 
ability through channel n is finite (0 < T„ < 1), then 



there is a non-zero probability p„ = 2(1 — T„)T„ that 
exactly one electron out of a singlet pair is transmitted, 
and this contributes to spin fluctuations. 

Another interesting conclusion which we can draw from 
this work is that, when measuring non-commuting quan- 
tities with subsequent detectors, one should take into 
account the quantum dynamics of the detectors them- 
selves. This is because the decoherence time for the de- 
tectors, Tc can be larger than the average time between 
two subsequent counts, T5. Thus if one would put by 
hand the off-diagonal elements of the detectors' density 
matrix to zero after each count, which amounts to apply- 
ing the projection postulate, one would obtain the wrong 
result. We have shown that, in the system considered 
here, such a deviation from the naive application of the 
projection postulate is revealed by the fourth correlator 
of spin counts. 
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